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Abstract 

Motivated by a need to understand spin-momentum transport in CPP (current perpendicular 
to the plane) structures, a quantum field theoretical treatment of spin-spin interactions in ferro- 
magnets is presented. The sd interaction of the conduction electrons and the magnetic medium 
is treated non-perturbatively from first principles in real space. The localized magnetic moments 
also interact with each other through a Heisenberg exchange potential. To take into account 
correlation effects, a second quantized formulation is used. The semi-classical limit is taken by 
using a coherent-state path-integral technique which also allows us to go beyond a linear-response 
approach. We derive a set of coupled equations of motion for the nonuniform magnetization, the 
spin current and the two-point correlation functions of the magnetization. The rate of change of 
the magnetization is shown to obey a generalized Landau-Lifshitz equation that takes into account 
interaction with the conduction electrons. Within the relaxation time approximation it is shown 
that the polarization of the conduction electrons obeys a diffusion equation. The diffusion tensor, 
which has off-diagonal terms due to the sd exchange interaction, is now explicitly dependent on 
the magnetization of the medium. We also show that the magnetization fluctuations satisfy a 
diffusion- type equation. The derived equations are used in two illustrative examples. 
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I. INTRODUCTION 



Spin-spin interactions in ferromagnetic metals are of paramount importance in today's 
GMR recording heads. There is also currently great interest in the magnetic recording 
industry in using spin currents, instead of magnetic fields, to switch the magnetization in a 
writer device. In this case a polarized electronic current is needed, such that the net spin of 
the polarization becomes effectively another magnetic source which induces an interaction 
with the magnetic moments of the medium. One widely used approximation is to separate 
the degrees of freedom of the current from those of the local magnetic moment. This latter 
separation is not justified in conducting metals but it nevertheless produced good results 
in some cases. 1 This paper explores in detail the consequences for the spin accumulation 
problem in ferromagnets of assuming that the interaction between the conduction electrons 
and the local moments is of the sd exchange type. This interaction can give rise to what 
is now known as spin-momentum transfer (SMT) in magnetic multilayers. This latter 
mechanism has been predicted by Berger 1 and Slonczewski 2 and later verified experimentally 
by various groups.^ Other interaction mechanisms between the conduction electrons and 
the magnetization vector have been proposed since the Berger-Slonczewski work . 5 i 6 i 7 i 8 In 
previous work, the interaction of the polarized current with the magnetization has not 
been treated self-consistently. In fact the equations of motion were based on those of a 
similar system, that of a current interacting with magnetic impurities.— We believe that 
this approach is not suitable for transition metals. 

In this work, we start from a microscopic description of the conduction electrons and 
the ferromagnetic medium and then take the semi-classical limit to derive equations for 
macroscopic quantities of physical significance to experiment and other phenomenological 
approaches. Although the derivations are somewhat complex, one can go to the main results 
(e.g., Eqs. and |HH which are generalized Boltzmann-type equations) and see that the 
correct physics is contained in them. 

Our results include those of reference [D and in fact our treatment should provide a 
basis for the more phenomenological work. We do not assume that the magnetization is 
uniform as in previous work and we will focus mainly on the diffusive regime. We use many- 
body field theoretical methods to describe the system of magnetic moments plus conduction 
electrons. Even though only a single particle picture is needed, the methods we use permit 
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us to treat the magnetic part of the problem and the conduction electrons on the same formal 
level. This allows us to derive transport equations for the conduction electrons and the local 
magnetic moments and include relaxation effects without recourse to more phenomenological 
modeling. Exchange effects, which are important in transition metals, are also included 
self-consistently. Finite-temperature properties are naturally included through the use of 
a path-integral formulation of the problem.— Including spin-dependent interactions in a 
transport problem means that we have to deal with many indexes. Path-integrals are 
helpful with book-keeping and hence simplify the discussion considerably as compared to 
Ref. 0. Finally a path-integral representation helps in making consistent approximations to 
the derived transport equations. 

The paper is organized as follows. In Sec. II, we introduce the Hamiltonian of the inter- 
acting electron-magnetization system in second quantized form. In Sec. Ill, a functional Z, 
which we call the universal generating functional, is defined in terms of the density matrix 
of the system and external sources. This functional has a path-integral representation and 
will generate all possible correlation functions. For the local magnetic moments we adopt 
a coherent state representation which is most suitable for a semi-classical treatment of the 
medium such as presented here. In Sec. IV, we derive equations of motion for the mag- 
netization (a modified Landau-Lifshitz equation) and the spin accumulation. These latter 
equations are then used to derive equations for the correlation functions of the magneti- 
zation. We also show that non-uniform magnetization of the medium gives rise to a spin 
accumulation effect similar to that due to interfaces.— This is one of the main results of 
this work. In Sec. V, we show how to solve these equations perturbatively and derive the 
terms that give rise to the spin momentum transfer effect for a system with non-uniform 
magnetization. Finally, we derive a diffusion-type equation for the magnetization fluctua- 
tions where the diffusion coefficient is dependent on the exchange integral. In section VI, 
we show some applications of the derived spin accumulation equations in simple cases where 
the magnetization is nonuniform in the direction of flow of the current. In the last section, 
we summarize our results. 
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II. THE MODEL 



We start from a quantum field picture of the electronic current and the magnetic moments 
of the localized electrons in a thin slab of thickness comparable to the mean free path, Fig. 
[T] In the following we do not include explicitly an electric field, but we assume that it is 
part of the spin-independent part of the Hamiltonian. Its inclusion has been done in Ref. 
I 111 . For the magnetic medium, a Heisenberg Hamiltonian is assumed. We explicitly take 
into account only exchange effects among the Heisenberg spins. Other important effects 
such as the demagnetization field term can be added phenomenologically. The conduction 
electrons will be assumed to be in equilibrium with the magnetic medium since they relax 
much faster than the magnetic moments. This is the case in assisted spin-momentum 
transfer switching. In fact the magnetic moments can be regarded as in contact with a 
Fermi bath with spin degrees of freedom. In this case switching may be regarded as a 
dissipative effect accompanied by a shift in the magnon energies.— Our work is a natural 
generalization of the model used by Langreth and Wilkins^ to study spin resonance in dilute 
magnetic alloys. The conduction electron field \I/ S satisfies the usual anti-commutation 
relations, 



In this system, the electrons are treated as non-interacting, i.e. no Coulomb interaction, 
and hence the electron field can be expanded in terms of single particle wave-functions: 




(1) 




where s is a spin index (s = ±|)- The number density operator of the electrons is given 





(3) 



s.i > 
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where a s (al) is the annihilation (creation) operator for a particle of spin s. The total 
Hamiltonian of the system is composed of a conduction electron part, H s , a magnetic part, 
H d and an interaction part, H sd , 



H = H S + H d + H sd . (4) 

The conduction electron part in a ferromagnet, such as Fe, is due to 4s-type electrons. The 
magnetic part is however due to 3d-type electrons. For the free conduction electrons we 
have 

H s = / 7^W f (r) • W (r) dr (5) 
where ^ is the two-component Fermi field, 




The magnetic medium is microscopically a lattice with a spin vector Sj at each lattice point 
i. Since we are interested in the continuum limit of this model, we can define a macroscopic 
spin vector for the medium 

N 

S (r) = £ S t 5 (r - r,) . (7) 

i=i 

The magnetization vector is then simply given by the volume average of the global spin 



M(x)=^ B ^, (8) 

where g is the spectroscopic splitting factor and /x_b is the Bohr magneton of the electron. 
The spin vector has an SU (2) representation and consequently so does the magnetization 
vector, e.g. 

[M x (x) , My (x')] = 2t 9f i B M z (x) 5 (x - x') . (9) 

We employ units such that g/iB = h— 1. Hence the magnetization will have opposite sign 
to that usually defined in the literature. The Hamiltonian of the spins is taken to be of the 
Heisenberg type, 

H d = -\j2 J ^- r ^ s ^^- B -J2 s - ( 10 ) 

ij i 
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We take account of exchange only and assume the spins to be in an external field B. In the 
continuum limit, we can write the exchange term in an integral form, 

i N if 

- Si ■ 8, J (r* - tj) = - / dvdv 1 S m to J (r - r') S„ (r') . (11) 

i,j=l J M=l,2,3 

Contributions from the magnetostatic field can also be included. For the typical thin 
structures we are interested in, the magnetization of the medium is usually held in-plane 
by the magnetostatic field. The energy term associated with this field is usually non-local 
but for a thin film it can be well approximated by a local term proportional to the normal 
component of the magnetization, 

H demag = l -K y rfx [S (x) • n] 2 , (12) 

where n is a normal unit vector to the surface and K is a constant depending on the shape 
of the slab. This latter term will lower the symmetry of the problem even further. For 
simplicity in the following we leave out this term and discuss it elsewhere.— 

The interaction between the electrons and the localized spins is taken to be of the s-d 
type, of the form 

H sd = ~ J dx (*t ( x ) ^ ( x )) . S (x) (13) 

where A is a coupling constant of the order of 0.1 eV and ~~a is a vector whose components 
are the Pauli matrices, 

[a i ,a j } = 2ie ijk a k . (14) 

e yfc is the antisymmetric unit tensor. The full Hamiltonian is then the sum of all the above 
terms, 



H=J2[ rfr^(r){-i-p 2 (r)+t/(r)-ia.BU a (r) 

a=l,2*' ^ J 

2 3 

a,/3=l i=l 

- ~ J dvdv' { S l (r) J (r - r') St (r') + 2^ (r) 5 (r - r') S t (r') } 



(15) 



i=l 
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where we have included a spin-independent external potential U (r) for the conduction elec- 
trons. The external magnetic field B can be taken to be time-dependent and/or spatially 
dependent in the following treatment. 



III. THE UNIVERSAL GENERATING FUNCTIONAL 



In this section, we introduce what we call a universal generating functional from which 
we derive equations of motion for the magnetization and the spin current. This functional 
is defined in terms of the density matrix p of the current-medium system and can be used 
to study equilibrium as well as non-equilibrium properties.— We are at present interested 
only in the near-equilibrium case. A general overview of the method appeared in Ref. 
Application of this method to spin systems is however almost absent. As far as we know, 
Refs. OandQ are the only work that apply non-equilibrium methods to spin systems. One 
of the other advantages of this method is that it allows equal treatment of thermal effects 
and non-thermal effects. The extension beyond a linear-response approach is also easily 
achieved, at least formally. This becomes essential when we are interested in questions that 
involve switching of the magnetization.— We will give a general outline of the method as we 
apply it to the particular s-d exchange system since the steps used in the derivation of the 
equations of motion are hard to find in one place. The importance of these methods, which 
are hardly used in magnetism, is potentially very great. 

To motivate the structure of the functional we are about to introduce, we recall the 
structure of the density matrix elements. In quantum mechanics, we usually use the energy 
method to solve a problem. However, in transport problems, time-dependent methods 
are essential. Hence the evolution of the density matrix comes into play. One of the 
ways to calculate density matrix elements is through a path-integral representation.— We 
have recently used such a formalism to treat the problem of fluctuations and dissipation 
in coherent magnetization.—*^ The density matrix elements at a given time t are usually 
written in terms of those at an earlier time using propagators running from past to present 
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and from present to past, 

(*\p(t)\&) = <&'(t)} (16) 

S)$i2)$ 2 ($ (t) |$1 (t )) ($1 ^O) |P| $2 (to)) ($2 (to) |$' (t)) 

D$!2)$ 2 ($! |p (to) I $2) <3 (+) ($ (t) , $1 (t )) <g H ($2 (to) , (t))(17) 



The propagators an d l2»(~) are then written in terms of path-integrals. This is equivalent 
to using the Schroedinger equation to solve for the evolution of the density matrix. In a 
transport problem, we instead introduce a functional of the density matrix. This functional 
is then made to depend on new virtual sources rji, rjl, r] 2 , r}\ , Ji and J 2 . These sources are 
coupled to the conduction electrons' field and the magnetic moments of the medium which 
will enable us to generate all kinds of correlation functions and their time-evolution. The 
functional is then given in terms of a trace formula, 



Z [ 771, 37*, r]2, 772, Ji, J2, P 



(18) 



Tr I p (t ) (V^exp -i J dx ( 77* (x) • * ( x ) + ( x ) - r, 2 (x) + J 2 (x) - S (x) ) ) 

| dx (r ?1 *(x)-*( X ) + * + (x)-r /l (x) + J 1 (x)-S(x)) ) } 



x ( Texp 



where T is the time-ordering operator. It orders operators with the earliest time argument 
to the right. T _1 is the inverse of T. The external sources 771 and r] 2 are two-component 
classical (i.e. Grassmann) sources where rji (772) and r)l (rj^) are treated as independent. 
The operators are all written in the Heisenberg representation. The need for both time- 
ordering operators is clearly seen through the Feynman- Vernon Influence formalism, Eq. 
(|16|) . The coefficients of the expansion of the functional Z in terms of its arguments give 
all possible correlation functions of the system. A few of the lowest order correlations are 
stated below. They can be easily verified by differentiating the functional Z with respect to 
the virtual sources. For example, to get the expectation values of the conduction electrons' 
field, we differentiate Z with respect to 771 or 772, 



Z5r] 



6 -^—\ v=r=J=0 = i(T¥ A (x)> = i<¥u(x)> , (19) 



1A 
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Similarly to get expectation values for the magnetization, we differentiate with respect to 
the external sources Ji and J 2, 

' ^ = i(r5i(x)) = i(5 u (x)) (21) 

r]=r]*=3=0 



Z 5 Jij (x 
1 <JZ 



= -z(T- 1 5 i (x)) = -2(^(x)). (22) 

JJ=JJ*=J=0 



Z £ J 2i (x 

Higher order correlations can be obtained simply via higher order differentiations. 

The Hilbert space for the conduction electrons and the local magnetic moments of the 
medium is the product of the corresponding Hilbert spaces, 

|*,n) = |$) ® (23) 

where |$) is a many-body fermionic state representing the conduction electrons and \fl) is 
a magnetic moment state. The magnetic moment states will be represented in terms of 
spin-coherent states (SCS)^ (and references therein). Since the operators are initially taken 
to be in the Heisenberg picture, then in the presence of the additional external sources, rji, 
r] 2 , and J, the states are no longer time independent. Now we write the functional formula 
in terms of a path-integral. Since we are not interested in the transient behavior of the 
interaction between the current and the magnetic moments, we assume that the external 
electric field was turned on a long time ago and we will eventually set to — —00. Reference 
I12I treats the case where to is kept finite in a finite size thin film. Moreover, we assume that 
the density matrix is separable initially, i.e., 

p(t= -00) = p s (t = -00) ® p d (t = -00) , (24) 

where p s is the density matrix of the conduction s-electrons and pd is that of the local 
magnetic moments. Now let Oj) be an initial over-complete set of states for the operators 

^ (r, t ) <8> 1 and 1 (g> S (r, to). Similarly, we let |$ c , f2 c ) be an over-complete set of states 
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for the operators ^ (r, t c ) <8> 1 and 1 (g> S (r, t c ) at the time t c . At each intermediate time, 
we define similar states. Then the functional Z can be written as follows 



Z[»7i, r)i, V2, vh J i) J 2, P 



T 1 exp 



y rfx (r ?2 *(x)-*(x) + * + (x)-7 ?2 (x)+J 2 (x)-S(x)) 



x($ c ,O c 



Texp 



y dx ( H (x) ■ * (x) + (x) • m (x) + J x (x) ■ S (x) ) 



(25) 



Hence we can now formally write the functional Z as a time-ordered path-integral around 
a closed path in time starting at t = t Q , passing through t = t c and then going back to 
t = to (see figure 2). This functional then has a path-integral representation similar to the 
equilibrium cas o 10 ' 21 , 

Z[r]*, r), J, p] = J 3**3*3m exp{iA[^*, m, 77*, 7?, J] } (26) 
x (* 2 ,m 2 I p I *i,mi) 



where we have used the following notation for the now classical tensor fields m and 



m = (m + , m_) (27) 
* = (* + , $_) (28) 

where + and — stand for the component that is propagating forward and backward in time, 
respectively. The field m is therefore a 2 x 3 tensor, while $ is a 2 x 2 tensor. Similarly, 
we write the source terms in terms of tensors. Using the following notation for the external 
sources, 

V = (Vi, m) , (29) 

J = (Jl, J a ) . 
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r] becomes a 2 x 2 tensor and J a 2 x 3 tensor. This notation greatly simplifies the 
manipulation of the path-integral. 

The action A is given as the difference of two actions; one due to the fields propagating 
forward in time and the other due to fields propagating backward in time, 

A [**, m, rf, 7j, J] = A [**, mi, 77*, Vl , J x ] - A [**., * 2 , m 2 , 77*, r/ 2 , J 2 ] . (30) 

Both terms on the right are obtained in the usual way. The electron contribution is standard. 
The magnetic moment contribution can be obtained in the same way, but it involves a 
geometrical part coming from the SU (2) symmetry. Hence the forward part of the action 
is given by 



A [**., *i, mi, t£, 7]t, Ji] = A W z [mi]- J dxH d (mi ( x) ) (31) 

+ J dx jz a (x) ^*i a (x) - F s+Sd (*} a , *i a , mi) | 

+ y {rilaW *i « (x) + *J a (x) 771 Q (x) + J x (x) • m x (x)} , 

where a summation over a, the spin index, is implied. The Awz is the Wess-Zumino part 
of the action A. Because of the boundary conditions on the spin fields at the left ends of 
the time path at t — — oo (KMS-type conditions), this WZ-term has the same form as in 
the equilibrium case where the path of integration is along the imaginary-time branch from 
t = to t = -ip^Z. 

A W z = [ dr [ dt m(t, r) • [d t m(t, r) x d T m(t, r)] . (32) 
Jo Jc 

The vector map m(£, r) is a parametrization of the surface enclosed by the trajectory of the 
magnetization, Fig. |2] 

m(t,0) = xa x {t),te Ci (33) 

= m 2 (t), t G C 2 , 
m(t, 1) = m , 
m(— oo + z0 + ,r) = m(— oo + iCT, r). 

(34) 

m is a distinguished vector and is usually taken along the quantization axis. 
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This WZ-term is topological in origin and can be considered as a constraint condition 
on the configuration space of the magnetic moments and gives rise to magnetic monopole 
type potentials, 23 

Next we expand the initial density matrix elements in terms of the initial configurations of 
the conduction electron field and the magnetization field. Therefore we are led to define 
a new functional F which may describe any initial correlations between the conduction 
electrons and the local magnetic moments, 

(* 2 , m a | p | *i, mi) = exp { i F [* f , m] } . (35) 

Since we are assuming that the density matrix of the whole system is decoupled at t — to, 
then all cross terms in the expansion are zero. Keeping only terms up to second order, the 
expansion is 



F [m, * f ] = C (0) + J rfx e a/3 C^ (x) • (x) (36) 
+ l - J dxdy e Q W a (x) ■ C$ (x, y) ■ m p (y) 
+ J dxdy e Q V A *t (x) • Q 7 a (x, y) • * (y) . 

The tensor e is defined such that en = —622 = 1, and ei2 = £21 = 0. The functional 
coefficients C^ -*, C^ l \ C^ 2 \ and Q are as yet unknown. The notation used should be clear; 
for example the last term involves summations over the path index and the spin index, 

e Q V A *t . Q 7A . 9p = e^e^t a Q^'^ (37) 
where the upper indexes on Q are for spin up and spin down. 

Inserting this expansion back in Eq. (|26[). we end up with the following expression for the 
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action A, 



X[^*,m,77 + ,77,J,Q,C] =e a ^A wz {mp}- j dx H d {rap ( x) ) J 
/ rfx (z *J (x) (x) - tf s+sd , ^, m/3 



J fl fx) • m« (V 



(38) 



| rfx [r,l (x) ^ (x) + *t (x) ^ (x) + 

+ l e -7 e /3A y rfxdy mo (x) . C ^ A (X; y) . m/3 (y) 

+e Q V A y rfxrfy *t (x) . (x? y) . ^ (y) ? 

where we have made an obvious redefinition of the coefficients. The functional Z, now 
becomes of the standard form 21 

Z[^^,J,Q,C] = (39) 
y 3**3*3m exp { i .4 [* f , m, r/ 1 ", rj, J, Q, C] } . 

The integral notation emphasizes that the path in time is closed, Fig. [21 Therefore we now 
can apply the usual field theoretical methods to extract the equations of motion from this 
functional. From the correlation functions, it is clear that the functional 

W [77+, 77, J, Q, C] = -i InZ [ V \ 77, J, Q, C] (40) 

is the generator that we need to derive the irreducible Green's functions of the system. To 
get the average value of the conduction electron field or the magnetization, we differentiate 
with respect to the coefficients in the linear terms. For the conduction electrons, we have 



Srfa s (x) 

and for the magnetization, we get 

5W 



^<**.(*)>, (41) 



e a V^(x)). (42) 



SJai (x) 

The average of the conduction electrons field, a Fermi-type field, is set to zero while we set 
the average of the magnetization to be 

M ai (x) = (m ai (x)). (43) 
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Given the above definitions, Eqs. (|41l42j) . the two-point correlation terms are easily obtained, 



5Z 



5\n 



• e^e*<*t.(x) (y)>, 



Z5Q-'(x,y) 6Qi°'(x,y 



5W 



^ A (x,y) 



(44) 
(45) 



™ t aa 'eW'{*l s (x)*P'Ay)}, (46) 



% s (x) 8r)l s , (y) 

where s, s' are for spin up and spin down and i, j are for the spin field components. The 
indices a, (3 denote the branch of time in fig. 2. Mixed correlation functions can be 
obtained in the same way. 

Clearly, solving for the two-point propagators is the least we can do to have a meaningful 
solution that includes relaxation effects. Knowing these propagators amounts to knowing 
the particle density, the spin density, the current density, and the scattering amplitudes, 
among others. Since we assume that the conduction electrons' field has no mean value, its 
two-point propagator is then explicitly given by time-ordered products, 



% -, (x > y) = <T- 1 (*.(x)*+(y))>, (47) 

gg (x,y) = <*. (x) *J (y)>, 

G{( (x ) y) = <T(*.(x)*+(y))> > 

Ql( (x,y) = -<¥+ (y) * s ( X )>. 

From the above expressions, it is clear that the function Q21 is the "less than" Green's 
function and Q12 is the "greater than" Green's function. Q\\ is the Feynman propagator, 
while Q22 is the Dyson propagator.— These Green's functions are not all independent. From 
their definitions, we can see that 

Q{( (x, y) + g£ (x, y) = (x, y) + g s 2 ( (x, y) . (48) 

The Green's function Q12 is of special interest since it is related to the average of the density 
operator of the conduction electrons. The two-point functions for the magnetization are 
similarly given by 
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•M|(x,y) 
-MS(x,y) 



(T- 1 (^(x))(5,(y)))-(^(x))(^.(y)), 
<^(y)^(x)>-(^(x)><^(y)>, 
(r(^(x)^(y)))-(^(x))(^(y)), 
(^(x)^-( y ))-(^(x))(^( y )). 



(49) 



These Green's functions are easily related to the retarded and advanced Green's functions. 
Since we are considering a situation which is not far from equilibrium, we will follow closely 



the treatment in Ref. 



25l Therefore, as in the equilibrium case, we relate the "less than" 



functions to the distribution function of electrons and spin in the semi-classical limit. 



IV. THE VARIATIONAL PRINCIPLE: EFFECTIVE ACTION METHOD 

Since the functions J, Q, C are not bound to a simple physical interpretation, we make 
the following Legendre transformation, 



r 



M Q4 (x),^'(x,y),^(x,y) 



a/3 



w 



t r\ ss r^- 7 

'Jai, ^a/3 



- / dx J ai U)M a 



(50) 



dxdy Q% (x, y) <?£ (y, x) 



--J dxdy (x, y) (M& (y , x) + M/s j (y) M a , (x)) 

We end up with a functional T that is expressed solely in terms of magnetization and 
correlation functions of the current and the localized spins. The equations of motion are 
then found by differentiating T with respect to its arguments, 



5T 



<5M Q4 (x) 
ST 

(x, y) 
ST 



"Jai (x) 



-Q s *(y,x) 



dyC%{^y)M M {y) 



-A (y,x) 



(51) 
(52) 
(53) 
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Using the standard tools of field theory2£, we solve for J, Q, C in terms of M, Q and A4. 
A discussion of Wick's theorem is beyond the scope of this paper. Omitting terms of order 
A 4 and higher, we have the approximate effective action for the conduction electrons and 
the localized magnetic moments, 



T\M,g,M] = A [M] + -lndet AT 1 

S 2 An 



(54) 



+ - / dxdy 



+ / dxdy 



5M ai (x.)M Pj ( y y ' 
S 2 Ao 



i In det Q 1 



6W a ,(x) 5^ S ' (y) 



+ 0(A 4 ). 

The functional Ao is the functional A with all the source terms set to zero. The tensor g 



is equal to 1 if i — j — k — 1 and equal to -1 if i — j — k — 2 and zero, otherwise. The last 
term, which is clearly valid for large magnitude of S, has a simple interpretation in terms of 
Feynman diagrams, Fig. El 

The equations of motion for M, Q, Ai, are obtained by minimizing T and setting the 
external sources to zero with the appropriate initial conditions. Within the above stated 
approximations, the magnetization of the medium obeys the following equation of motion 



(55) 



e a0 e ilk (x) g tMpk (x) + e «/3 5a(] ^ (x) 

+ \ J e a/3 V 2 M^ (x) + | e a ' a a\, s g s J a (x, x a ) = 0. 

Here we have taken the long-wavelength limit to get the familiar exchange term through a 
coarse-graining procedure where each cell is taken to have a maximum spin of S. The last 
term on the left is clearly the interaction with the conduction electrons' magnetic moments 
to all orders in A. The equation of motion for the conduction electrons is 



(id ty - e (y)) 5 S , S + ^I'sB 1 + ^a ss >M ai (y) 



A 



(56) 



i an f.. o & <,u 



2 2 



G3$&*)G% (x,y 



a"/3" 



-i5 s 3 



7/ 3 i z -yj 
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where G (y) is the spin-independent energy of the conduction electron. The term of first 
order in A describes the full exchange interaction between the magnetic moments of the 
localized electrons and those of the current. The structure of this equation is familiar from 
the theory of correlation functions due to Coulomb interactions^ 27 There, the propagator 
M. (x, y) is replaced by the Hartree propagator, 
should follow by analogy with the treatment in Ref. 
of motion for the magnetic correlation functions, 



?herefore the solution of this equation 
1 

271 The final equation is the equation 



e aP e ijk {dtMpk (x) _ ^ (x) Qt} M wj (Z; x) + e ap I dy J ( x - y) M.% (Z, y) 



,k'j 



aj3 



a: 



dy 



iSli (x - 



The integrals are all four-dimensional and hence we have defined J (x — y) = 
J (x — y J 8{t x — t y ). This set of coupled equations, Eqs. (fo31 - \ 57)) . is the main re- 
sult of this work. It is important to observe that up to this point, the propagators Ai and 
Q are the true propagators of the theory. Hence the above equations are non-perturbative 
in nature. In this section and the next, we study the structure of these equations and make 
contact with previous work. Since we are interested in how the magnetic moment of the 
current influences that of the medium (or vice- versa), we define 



QJt^(x,y) = -o^Qil (x,y) , (58) 

to be the conduction electron spin propagator. The spin "charge" of the current is easily 
seen to follow from Wl a /3 by setting a — 1, (3 — 2 and letting y — > x + , 

97l 4 (x) = ^ s ^(x,x). (59) 

However, we will find it useful to go to the center of mass and relative coordinates (Wigner 
coordinate system) to make contact with the classical quantities, 

X=~(x>+y% (60) 

T = — (t X + ty) , 

xa = x" - "y, 

t A t x ty . 
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In this new coordinate system, the magnetic moment of the conduction electrons 97T(x) 
becomes a function of the macroscopic variables X and T only, 

OT(x,y) = J^-J ^exp[^ A -2p-x A ]9jr(X,T;cu,p). (61) 

To get the equation of motion for 971, the spin charge of the conduction electrons, we first 
multiply Eq. ()56)1 from the L.H.S. by a l ss „ and sum over the spin degrees of freedom. We 
end up with an equation for the polarized current propagator, 

ie a/3 

e af3 5 SS > [id ty - G J Wl ja (y , z) H — — [B + AM (y)] x 97t 7Q (y , z) (62) 
+AV° V '/ /3 ' /3 "%% / rfx[^-(y,x)>l^(x,y)^ v < s (z,x)] =0. 
The last term on the left provides for the relaxation of the spin moment of the conduction 
electrons. 

To derive the equation for the polarization of the current, in the following we use the 
relaxation time approximation and replace the last term in Eq. (J62|) , the collision integral, 
by a local term. The classical polarization of the current 97t c is found by first assuming 
that the I— th component , 971' (x), has the following form 



Wl l (T, X, p, oj) = 6 (u - p 2 /2m - V (T, X)) (T, X, p) , (63) 

where we have set G (p) = p 2 /2m + V(T,X). Then by averaging over the fast degrees of 
freedom, we have by definition 

m[(T,X) = v J ^-^Wl l (T,X,p,u;), (64) 

where v is the volume of the system. The spin current $ is defined in the usual way. However 
here it has a tensorial character because of the vector character of the spin charge, 

v kl frr> v\ V f myk frp v \ I 



zkl = _ , _^ ypp x ^ (65) 

m J [Ztt) 6 

The equation of motion for 9Jt c is found by first going to the center of mass coordinates and 
using the quasi-particle approximation. 24 We find for each k— component 
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d T + ■ d x Wl k (T, X, p) + e klp \B l + XM l (T, X)l Wl p (T, X, p) 
m J 

aji fc (t, x, P ) - wi k q m k (t, x, P ) - m k 



(66) 




where t&(t p ), is the relaxation time for spin flip (momentum) scattering processes . 2 ^ 2 ? By 
definition the average of the last term over the momentum is zero. This way of writing 
the collision term is valid only in the absence of spin-momentum coupling terms such as 
L • S-coupling. 

It is worthwhile to pause here and consider the content of this equation. The first term 
on the left hand side is the total time derivative, with independent variables (T, X, p) and 
dp/dt = (we consider a non-zero electric field elsewhere^ 2 ) . The other term on the left 
hand side is the torque due to the local moments and the right hand side is an approximate 
expression for the collision operator. 

In the absence of gradients, the polarization of the conduction electrons is along the 
local effective field and there will be no spin currents. This will not be the case if there 
is an electric field present. The spin accumulation effects studied here are solely due to 
non-homogeneous magnetization of the medium. It is important to observe that this effect 
is present even in the absence of a current as is the case for the spin accumulation effects 
due to interfaces put forward by Berger.— In fact paramagnetic-ferromagnetic interfaces are 
naturally included in our treatment since we are dealing with nonhomogeneous magnetiza- 
tion. An interface simply corresponds to an abrupt change of the magnetization from a 
non-zero value to a zero value. These particular effects will be treated below. 

Next we multiply Eq. ()66|) by the velocity and then average over it. To obtain the usual 
Fick's law for spin diffusion, we assume that the momentum relaxation time is small and 
hence the R.H.S of Eq. ()66j) is larger than the effects due to the local magnetization M. In 
this case we have 



where the diffusion constants are given in terms of an average Fermi velocity vp by 



3« (T;X) = _ D k dxiTl k ( TiX ), 



(67) 



nfe _ I„2 T e. 
_ 3 



(68) 



19 



where I/t^ = 1/7% + 1/v To get this result we made use of the following approximation 
for the velocities of the conduction electrons 

d Xl J ^v l vm k (T, X, p) ~ \v 2 F d X] m k c (T, X) . (69) 

Finally using Eqs. (j661 l6*7j) . we find that the classical magnetization of the conduction elec- 
trons obeys a diffusion equation for each one of its components, 



[&r - D k v 2 ] m k (t, x) = — (m k (t, x) - m k q (x)) - [(b + am) x m c } k . (70) 

This equation is however rotationally invariant and does not show the reduced symmetry 
of the ferromagnetic state. To get a more realistic equation we improve on Fick's law by 
keeping all terms in Eq. ()66)1 and treat exchange effects between the conduction electrons and 
the magnetization more carefully This amounts to taking into account the sd-exchange 
term in the electron propagators. For slow variations in time, we have now a modified Fick's 
law that takes into account the variation of local magnetization in space and in direction, 

Z kj (T, X) = -V kp d x m p c (T, X) , (71) 

where now the diffusion constant becomes a tensor. It is defined in terms of a matrix A 



£) fcp (X) = D p (A-^CX). 



(72) 



There is no summation over p in this equation. The matrix A depends locally on the effective 
magnetization field H, 



A(X) 



1 



TyH z 



H Z T X Hy 



7~zHy T Z H X 



(73) 



In our approximation, the effective local field is simply 

H = B + AM (x) . 



(74) 



Now in the steady state, the equation satisfied by the average magnetic moment DJl c becomes 
a generalized diffusion equation 



d Xl (x) d Xl m p (x)] = - (m k c - m k q ) + a [m x m c 



(75) 
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The tensor character of the diffusion term in Eq. (|75|) is not due to anisotropic transport 
- the flux in the j— direction in Eq. (|7Tjl is due to a gradient with respect to Xj. Rather, 
the p-component of DJl c is rotated into the k — th direction by the effective field H, while 
transport takes place in the direction of the gradient. The diffusion tensor, Eq. (J72|) . has 
striking similarities to the diffusion tensor of charged species in a plasma.— If we restrict 
ourselves to the case where the local effective field is constant and along the z— axis only, 
then the transverse diffusion coefficients are similar to those found by Hirslj^ and Kaplan^ 
using very different methods from the one presented here. Their work showed that in the 
direction perpendicular to the effective field, diffusion of polarization of the electron gas is 
much slower than along the field. The off-diagonal terms have their origin in the strong 
exchange interaction among the conduction electrons which can not be treated perturbatively 
for a transition metal. In a magnetic metal such as Ni, the off-diagonal terms can be two 
orders of magnitude larger than the diagonal ones.— 



V. SPIN-MOMENTUM TRANSFER: A SELF-CONSISTENT TREATMENT 



To solve the above equations of motion, we retain a subset of the terms arising in the 
full spin propagator. This approximation is essentially similar to the random phase ap- 
proximation in the calculation of the ground state energy of an interacting electron gas.— 
The zero order propagator is taken to be that of the electrons in the external B-field and 
the localized spins interacting through the exchange interaction in the presence of the B 
field. Since the magnetic moments constitute a many-body problem, a full solution is not 
possible in general. Hence an explicit solution to the problem requires first a calculation of 
the background magnetization. The magnetization satisfies a generalized Landau-Lifshitz 
equation which follows from Eq. (|55|). First we observe that when the external sources are 
turned off, we have 

Mi(x) =M 2 (x) = M(x), (76) 

where M is the average, i.e. classical, magnetization. Equation (|55|) is a system of two 
equations for Mi and M2, the magnetization vector along the paths G\ and C2 respectively, 
Fig. |21 It is the averaging of these two equations that gives the equation of motion for the 
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average magnetization, 
d t M (x) = M (x) x 
where 



l -J V 2 M (x) +B + ^ oV.± (pit (x, x + ) +gg (x, X") 



(x,x+) = < ( X , y ; 



(77) 



(78) 



The last term is simply the spin of the current. Recalling that at equal times and equal 
positions, all different Green's functions are related, the equation of motion for M simply 
becomes 



d t M (x) = M (x) x 



■J V 2 M (x) + B + A 9Jt (x) 



(79) 



The last term gives rise to dissipation and a contribution to the precession for magnetic 
multilayers.— As we will see below, this term becomes J-dependent in the non-uniform case. 
This latter equation, Eq. (|79|l . is the equivalent of the Landau-Lifshitz equation (LL) for the 
magnetization in the presence of a current. This form is still valid even in the presence of 
an electric field. The solution of Eq. (J56|) can be represented in terms of Feynman diagrams, 
Fig. 0] First, we define the propagator of a non-interacting electron in an external magnetic 
field B and zero electric field, 

0$ a " (*, y) = { K- E a ) + la*,!? } 1 SS (x, y) • (80) 
and expand Q in powers of A using Eq. ([56)1 . Keeping only terms up to order A 2 , we have 
0(0) -\ a g(o) M . a g(o) + lAV^M- ff ^M-ae (0) + X 2 G^a ■ M-aQ^aQ^ ... (81) 

This is a matrix equation and hence integrations over time, space and spin degrees of freedom 
are implicit in the above notation. Recently Mills calculated the damping contribution to 
order A .^ 12 One of his conclusions is that this contribution is dependent on the symmetry 
of the system in this case. This follows from the fact that the spin propagator of the 
conduction electrons in a ferromagnet is not 0(3)-invariant, to first order in A and higher 
since it depends explicitly on M. 

Equation (|57|) gives the dependence of correlations on the exchange interaction and on 
the s-d interaction between the current and the medium. Our assumption is that exchange 
interactions are much stronger than the spin-spin interaction. Hence to lowest order, we 
neglect the latter in the equation for the fluctuations. To understand the meaning of such 
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an equation, we study the case with strong exchange interactions, i.e., we take the average 
magnetization to be a constant and assume the external B field to be small. In this case 
Eq. (J57J) becomes 

d ty M* k (y, z)- [M ■ d ty M* k (y, z)] M = iff* (nxM)+Jdx [J (x - y) M% (x, z) x M] , 

(82) 

where for each k, the unit vector n has components n % = 5 lk . The notation A4 xk is for a 
vector with components A4 tk , i — 1, 2, 3. Now if we average over the variable z, we get an 
equation that gives the time variation of the fluctuations of the magnetization around M. 
These fluctuations will in turn cause fluctuations in the current through the last term in Eq. 
(|81|). Next we show how this latter equation gives rise to a Boltzmann-type equation for 
the magnetization fluctuations Ai d . First we expand A4 xk (x, z) around the position y , 



M xk (x, z) = M xfc (y , z) + d x M xk (x, z) | x=y Ax+ ^d^M xk (x, z) 



AxAx+... (83) 

x=y 

where Ax = 3? — y . If we put this back in Eq. (|82|). we get a diffusion- type equation for 
all components of the magnetization fluctuations 

dt y M x a k (y, z) - [M ■ d ty Ml k (y, z)] M + l -J 2 (y) M x VjA<J* (y, z) = ie^n x M (84) 

- J (v)MxM$(y,z), 

where Jo (y) and J 2 (y) are the zeroth and second moment of the exchange coupling, 

J (y) = y rfxJ(x-y), (85) 

J2(y) = lj dxJ(x-y)Aa; 2 . (86) 

These integrals converge since the exchange coupling is short ranged. The first moment van- 
ishes since we are assuming isotropic exchange coupling. Hence treatment of the coupling of 
the conduction electrons to the ambient magnetization at low temperatures or temperatures 
close to T c must include Eq. (|84jl to account for the fluctuations of the magnetization. 

VI. APPLICATIONS 

In this section, we mainly show how this formalism can be applied to multilayers. In Ref. 
fllL we showed how our results extend those of Zhang, Levy and Fert 7 by taking into account 
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the indirect exchange effect of the magnetization on the conduction electrons. This is an 
important effect in transition metals and can not be treated by a Born approximation. In 
the following we study two types of structures with nonhomogeneous magnetization. First 
we examine CPP-type structures with very thin paramagnetic spacers and no interfacial 
scattering. Second we consider structures which are topologically equivalent to a torus. 
These examples clearly illustrate the origin of spin accumulation to be directly related to 
inhomogeneities in the magnetization. It is also obvious from these examples that domain 
walls are another physical example where the results presented here can be applied. The 
interface will not be represented by a step function in the examples below and will instead 
take the shape shown in Fig. which plots the mean field a(x) due to the magnetization. 

For a local magnetization which is a function only of distance x in the direction of the 
current, M = M (x) z, the spin accumulation obeys the simplified equations 



d 2 m x (x) d 2 m y (x) 

dx z dx z 

2 Ply da (x) dm x (x) | 2 Dxy \ ^ — dm " y — - — Q — my — (87) 

Pa (x) dx dx \ xx P J dx dx r s f T s f 

d 2 m x (x) d 2 m y (x) 

xy dx 2 + yy dx 2 
n n Pl y \ da (x) dm x (x) P 2 xy da (x) dm y (x) _ m y (x) t a(x)m x (x) 



P J dx dx Pa (x) dx dx r s f r s / 

^d 2 m z (x) m z . . 

D ^ A = m 

where 

a (x) = tXM (x) . (90) 

The coefficients P xx , P yy , P xy and P yx are functions of the local magnetization, the scat- 
tering rates and the exchange constant: 

P 

Dp Pxx Pyy " . ' ~ , 7~oi (^^) 

1 + (rAM(x)) 

Pyx = 'Pxy = 1 + (rAM(x)) 2- ( 92 ) 
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These equations will be solved for different configurations of the magnetization M. We adopt 
the following parameters for our calculations: The spin diffusion length l sc u = ^Dr s f = 
100 ran, D = 10~ 3 m 2 /sec, D yx = 100 D xx , and A = 0.1 eV. It should be noted that in these 
equations, the torque term has the opposite sign to that which appears in e.g., Zhang et 
al.-i since we have taken the electron charge to be positive in our definitions of the magnetic 
moments. 

First we consider a configuration with in-plane magnetization. The magnetization is 
assumed to vary with position along the direction of the current. The spacer has practically 
zero thickness, which is a reasonable approximation for most GMR devices. Fig. shows 
a(x) for a typical interface. We do not explicitly include a non-magnetic spacer but we 
set the magnetization to zero at the center. At the ends it is parallel to the local z-axis, 
the local direction of the equilibrium magnetization. The transverse components of the 
spin accumulation are set to zero at the outer ends. The respective z-axes make a non- 
zero angle in order for the spin accumulation to be non-zero. First we demonstrate the 
effect of inhomogeneities on the spin accumulation; we keep the relative angle between the 
magnetizations the same but we vary the size of the 'domain wall', the transition region of 
the local mean field. It is seen from Fig. IHland Fig. Qthat the larger the inhomogeneities the 
larger is the spin accumulation. This effect is independent of the relative orientations and 
was not predicted before. Therefore spin accumulation can be enhanced by having a layer 
with constant direction magnetization but with spatial inhomogeneities. Such a structure 
can be achieved by e.g. having a temperature gradient across the slab or controlled doping 
that changes the magnetic saturation along the direction of the current. Figures El and El 
show the variation of the z-component of the spin accumulation with respect to the relative 
angle of the magnetizations and size of the sample. Figures El and ^2 ar e for the x and 
y components for the smaller sample. Here, the spin accumulation is largest for the case 
where the two magnetizations are orthogonal to each other and the size of the sample is 
smallest. In all these results, the equilibrium spin accumulation is normalized to -1 on the 
left hand side and normalized to +1 on the right hand side. The components of the spin 
accumulation are taken with respect to a global frame, that of the layer on the left. 

As our second example, we choose a 'ring' structure. This can be part of a long solenoid 
with a square cross section. Hence we now solve our equations with periodic boundary 



25 



conditions. In each side of the square cross section, the profile of the magnetization within 
a period is shown in fig. The relative angle of the magnetization between neighboring 
sides is 90 degrees. We study the spin accumulation as a function of the size of the sides. 
In figs. El and E3 we plot the three components of the total spin accumulation. The 
equilibrium spin accumulation is taken to be normalized to one. The solutions show the 
expected behavior. The spin accumulation tries to reach its equilibrium value near the 
middle of each side. The spin accumulation is largest when the length of each side is 
smallest, as expected. This geometry shows how spin accumulation can be transported 
over large distances and also modulated by controlling the size of the cross section, similar 
to what happens in a regular transformer except here we are working with spin charge. 



VII. CONCLUSION 



We have introduced a many-body formalism based on path-integral techniques capable 
of handling a system of both local magnetic moments and conduction electrons in a self- 
consistent manner. Transport properties can be obtained through the calculation of the 
two-point functions of the current and the magnetization respectively. One of the impor- 
tant outcomes of this treatment is that we were able to derive a set of new equations that 
are needed when the magnetization of the medium is no longer homogeneous. First we 
showed that the polarization of the current is no longer homogeneous and satisfies a gen- 
eralized diffusion equation where the diffusion tensor is dependent on the direction of the 
magnetization. We have hence shown that exchange effects are important in a ferromagnet 
and need to be taken into account properly. The fluctuations of the magnetization were 
also shown to obey a diffusion-type equation which depends on the direction of the local 
magnetization. This latter equation is in addition to recovering a Boltzmann equation for 
the current which follows from Eq. (Jo^j) and a Landau-Lifshitz equation for the average 
of the magnetization, Eq. ([79)1 . We have also shown how the non-uniform magnetization 
affects both the conduction electrons and the spin-momentum transfer term. We gave a 
simple physical picture for our main results. We finally showed how our results can be ap- 
plied in various configurations. Our results show that spin accumulation can be enhanced 
by inhomogeneities at the interface. 

In this work we focused on spin accumulation and we did not deal with its effect on 
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the dynamics of the local magnetization. We showed how to recover the Landau-Lifshitz 
equation but we discussed the effect of the collisions on the classical magnetization only in 
qualiative terms. The use of the relaxation time approximation is another shortcoming of 
this work. Its improvement will complicate the treatment further. We believe these issues 
should be addressed in the future. 
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FIG. 1: Geometry of the magnetic sample. 



FIG. 2: Closed time path: branch 1 corresponds to forward propagation in time while branch 2 is 
that for backward propagation. 
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FIG. 5: Profile of the interface (or molecular field a(x)) used in the text. The current flows 
perpendicular to the interface. The nonmagnetic spacer is taken to have zero thickness. 
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L1 = 1 nm, theta = 45 




FIG. 6: The x-component of the spin accumulation as a function of the interface inhomogeneities. 
On the left, M = —Mqz and on the right M = M$z. Ida = 0.6nm corresponds to the sharper 
interface. 
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LI = 10 nm, theta = 45 




FIG. 7: The y and z components of the spin accumulation as a function of the interface inhomo- 
geneities for the same configuration of M as in Fig. HO 
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LI = 12 = 100 nm 




FIG. 8: The z-component of the spin accumulation as a function of the angle for L = 100 nm. The 
angles shown on the right are expressed in degrees. 
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LI = L2 = 10 nm 




FIG. 9: Spin accumulation as a function of the angle between the magnetization in both regions: 
z-component for L = 10 nm. The angles are in degrees as in fig. 8. 
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LI = L2 = 10 nm 




FIG. 10: x-component of the spin accumulation as a function of the angle between the magnetiza- 
tion in both regions for L = 10 nm. 
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LI = L2 = 10.0 nm 




FIG. 11: y-component of spin accumulation as a function of the angle between the magnetization 
in both regions for L = 10 nm. 
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FIG. 12: Profile of the local magnetization along the different sides of the square ring (or torus) in 
one period. The magnetization is in the yz-plane whereas the current is in the x direction. The 
relative angle of the magnetization between neighboring sides is 90 degrees. 
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LI = 10 nm 
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FIG. 14: Spin accumulation in a torus strucure with equal sides intersecting at right angles. The 
length of each side is L = 100 nm 
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LI = 10 nm 
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FIG. 15: Spin accumulation in a torus strucure with equal sides intersecting at right angles. The 
length of each side is L = 1000 nm 
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